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Abstract 



We compare the energies of the Laughhn Hquid and a charge density wave 

in a weak magnetic field for the upper Landau level filling factors i^n = 1/3 

and 1/5. The charge density wave period has been optimized and was found 

to be ~ 3Rc, where Rc is the cyclotron radius. We conclude that the optimal 

charge density wave is more energetically preferable than the Laughlin liquid 

for the Landau level numbers > 2 at z^at = 1/3 and for > 3 at i^at = 1/5. 

This implies that the 1/3 fractional quantum Hall effect cannot be observed 

for > 2, in agreement with the experiment. 
PACS numbers: 73.40.Hm 
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The fractional quantum Hall effect (FQHE) was first discovered at the lowest Landau 
level (LL) This remarkable phenomenon occuring at certain unique values of the filling 
factor V = 1/3, 1/5, . . . has been associated with the formation of a uniform incompressible 
quantum state, or the Laughlin liquid 0. The traditional alternative to the Laughlin liquid 
is a charge density wave (CDW), which does not exhibit the FQHE. The FQHE occurs 
because the Laughlin liquid is lower in energy than the optimal CDW, which at the lowest 
LL has the same spacial periodicity as the triangular Wigner crystal [^. 

Later, the FQHE was found at the first excited LL also. The theoretical work, driven by 
this discovery 0], addressed the question of what exactly ground state, liquid or crystalline, 
is formed at high LL's? The liquid state at the A^-th LL was defined as follows: 




Here a] is the inter-LL ladder operator, raising the i-th electron to the next LL, and 
is the Laughlin state at the lowest LL. As a crystalline state Ref. 0] continues using the 
Wigner crystal as well as at the lowest LL. However, a recent more elaborate investigation 
of the CDW state [^] shows that the optimum one should have a period of the order of Rc, 
the cyclotron radius. The formation of the CDW with this period, which is a characteristic 
spread of electron wave functions, enables the system to reach a lower value of the interaction 
energy. Below we compare the energies of this optimal CDW and the Laughlin liquid. We 
show that the crystalline ground state rules out the FQH states at high LL's. 

Our calculation is based on the following model 0. We explicitly consider only the 
electrons at the upper LL, which is assumed to be spin-polarized. All the other LL's are 
completely filled. The role of these lower LL's is reduced to the screening of the Coulomb 
interaction among the electrons at the upper LL. The screening is accounted for by means 
of the dielectric function [^]: 
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e(g) = l + v{q)U{q), 



m<N 
N<n 

nm \ 2 I ' 

where / is the magnetic length, k is the bare dielectric constant, v (q) = 27re^/Kg is the 
Coulomb potential, L^{x) is the Laguerre polynomial, and is the cyclotron frequency. 
This dielectic function tends to unity in the limits g — and g — oo and reaches its 
largest value of 1 + \/2Nrs at g ~ Rc^- Here = A/2//cpag, with k-p being the Fermi wave 
vector, and being the effective Bohr radius. This model correctly renders the low energy 
physics of the system in the limit <^ 1 and Nts ^ 1 [Ref. 0]. Moreover, the results 
obtained within the framework of this model remain correct to the leading order in even 
for Nts ^1. In the latter limit e (g) ^ 1, which is consistent with the fact that the LL 
mixing can be ignored completely. 

Let us now describe the CDW state at the upper LL. In the limit of a weak magnetic 
field (A^ ^ 1) a simple quasiclassical picture can be given. In this case electrons can be 
viewed as classical particles rotating in cyclotron orbits. The only constraint imposed by the 
Landau quantization is that the concentration of the centers of the cyclotron circles at any 
point does not exceed l/27r/^. One can fill in a disk with these centers at their maximum 
concentration (see Fig. p. We call this disk a bubble and the triangular crystal built out of 
such disks a bubble phase. It is shown in Ref. |0 that the optimum number of electrons in 
a bubble is 

M ~ 3z/^iV, (3) 

which corresponds to the separation ~ 3-Rc between nearest bubbles. The nonuniform 
distribution of guiding centers with such a period, chosen in accordance with the form 
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factor of the electron wave function, does not create significant variations of the charge 
density. Hence, the electrostatic (Hartree) energy of the system does not increase too much. 
However, the exchange interaction in this ferromagnetic state favors an increase of overlap 
among the wave functions and thus the most compact arrangement of the guiding centers. 

For the comparison of different trial states we use the cohesive energy, which is the energy 
per electron relative to the uniform state formed at high temperature. The cohesive energy 
of the CDW state has been shown to be -E^h^ = —rfioj^ oc —B [Ref. p], which is of the 
same order of magnitude as the exchange-enhanced spin splitting f^. 

As for the cohesive energy of the Laughlin liquid, arguments can be given that it is 
very close to the cohesive energy of the Wigner crystal. The lower bound for the latter was 
estimated in Ref. [|5[ to be i^^^j^ ~ —huj^/^/N oc —B^/'^. In the limit of small magnetic field 
the CDW should obviously be more energetically preferable, ruling out the FQH states. 

It is interesting to know at what LL the transition from liquid to the crystalline ground 
state occurs? To answer this question one cannot use the quasiclassical approach and the 
CDW should be defined more accurately. We do so in several steps. First, we introduce the 
wave function of one bubble consisting of M electrons at the lowest LL: 



Here Zj = Xj + iyj is the complex coordinate of the j-th electron. Second, we define the wave 
function of a bubble at the A^-th LL centered at point R by raising every electron onto this 



where bi is an intra-LL ladder operator. To finally obtain the wave function of the CDW we 
build an antisymmetric combination of the bubbles centered at the triangular lattice sites 





LL and shifting its position by the magnetic translation operator |^ 





^CDW 



Y.^ign{P)X{^>^^{P{T,)}. 



(6) 



4 



Here P are the permutations of electrons between bubbles. For the case M = 1 this trial 
state coinsides with the Maki-Zotos ansatz wave function for the Wigner crystal §]. It 
can be easily seen that ^^cdw is of the Fock type, and that the overlap between the wave 
functions of different bubbles is negligible. 

To show that ^E'cdw matches our earlier quasiclassical picture we introduce the guiding 
center density operator: 

Hr) = 2nPY.S{r-R^). (7) 

i 

The summation here is carried over the electrons at the considered LL and = r^ + j^z x Pj 
is the guiding center operator, with Pj being the canonical momentum of the i-th particle. 
The average of this operator has a physical meaning of the guiding center density |P] . It can 
be shown that for the state defined by Eq. 



(8) 



'M 

which is just the Fourier transform of a uniform disk with the radius W2M. The last 
equation in (H) follows from the asymptotic formula for the Laguerre polynomials [llO| . 

The cohesive energy of the CDW can be calculated in the same way as it has been done 
for the Wigner crystal [Q: 



^'^N q^O 



A 



(9) 



The summation in Eq. (|^) is carried over the reciprocal vectors of the triangular lattice. The 
Hartree-Fock interaction potential u-^p (q) is defined in the same way as in Ref. 0: 
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MhF = Mh (g) - Mex (?) 



(q) = 27r/2e (g) ^^^ ' (10) 



«ex(g) = 2vr/^/-^e'-''^«H(g')- 



(2vr)^ 

Using Eqs. (|^,§), the cohesive energy for any given u^- can be calculated numericaly. The 
result is, of course, different for different M (see Fig. Therefore, one has to find M 
corresponding to the lowest energy. The energies of the CDW optimized in this way are 
summarized in Tables | and 0. In Table | we present the results for the case Nrg <^ 1, 
when the LL mixing can be ignored completely, i.e., e{q) = 1. Table || shows the results for 
the practically important case = One can see that the optimum number of electrons 
per bubble is the same both with and without the screening and is in perfect agreement 
with Eq. (I). 

The above results have been tested by the self-consistent Hartree-Fock procedure, similar 
to that described in Ref . . Starting from the initial approximation given by wave function 
(|), this procedure finds the optimal set of {u (g)) for a given periodicity of the CDW. The 
obtained corrections are of the order of lO^^r^huJc and thus do not affect the significant digits 
displayed in Tables |I| and We associate the corrections with a slight nonorthogonality of 
the wave functions of different bubbles. 

Let us now discuss the Laughlin liquid at high LL's. The interaction energy per electron 
can be calculated using the density-density correlation function 

^"^"^ = W) ^''^ 

where PN{f) is the projection of the density operator onto the A^-th LL. This can be most 
effectively done in the Fourier space because (q) is very simply related to ho (g) (the 
correlation function for N = 0), obtained earlier by Monte-Carlo simulations Q]: 
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hN {q) = ho (g) 



2 



(12) 



The calculations are greatly simplified by virtue of the approximation formula for Hq given 



in Ref . 1TT| 



The cohesive energy per electron is then, according to the definition 

Ecoh = EcoT — E^^^ = EcoT + ^^ex (5' = 0) , (13) 

where E^^^ is the interaction energy per particle in the uniform uncorrelated electron liquid 
formed at high temperature. The results of the numerical evaluation of these energies are 
also listed in Tables | and |I|. 

Compare now the energies of the Laughlin liquid and the CDW. As one can see, at = 
and = 1 the Laughlin liquid is more energetically preferable. At large A^ however the 
CDW wins. The transition to the bubble state both with and without screening occurs at 
N = 2 for ujy = 1/3 and at A^ = 3 for z/jv = 1/5. The difference in the energies of these 
two states at A^ = 2 and 1/^ = 1/5 is very small. It should be noted in this connection that 
the trial state (|^), being compressible in nature, can be further improved by introducing the 
magnetophonon correlations ||12|. Another possibility arises when the system is subjected 



to an external impurity potential. The CDW can lower its energy by accomodating to this 
potential, while the incompressible uj^ = 1/5 liquid state cannot ||13[. Hence, the question 
of what exactly phase dominates at = 4| remains open. 

One can suggest the following interpretation of our results considering the zero-point 
vibrations of the lattice. Electrons at any A^ form the Wigner crystal when the upper LL 
filling is so small that the average distance between electrons is larger than the characteristic 
spread of electron wave functions Re- The amplitude of the zero-point vibrations of the 
crystal is mostly determined by the magnetic barrier since the interaction energy per electron 
is smaller than the cyclotron gap for the cases of interest. Hence this amplitude is ^ /. For 
the crystalline state to be stable this amplitude should be small compared to the lattice 



constant (the Lindemann criterion) . Now an important difference between the low and high 
LL's arises. At low LL's an increase of the electron density reduces the lattice constant. At 
some value of un it becomes of the order of / and the crystal melts into the Laughlin liquid. 
At high LL's, however, the lattice constant does not change as one increases the LL fiUing 
but remains of the order oi I. Hence CDW cannot be melted by quantum fluctuations 

at high LL's. 

In conclusion, we have compared the energies of the Laughlin liquid and the CDW with 
the optimized period (~ Rc) at the upper LL flUing factors vn = 1/3 and 1/5. We found 
that the 1/3 liquid state is unstable for iV > 2, while the 1/5 state loses to the CDW at 
N > 3. Our result imphes that the 1/3 fractional quantum Hall effect cannot be observed 
at flUing factors u > A. This conclusion is in agreement with the existing experimental data. 
The difference between the energies of the CDW and the Laughlin liquid at iV = 2 and 
i/jv = 1/5 is so small that more work is needed to distinguish them unambiguously. The 
authors are grateful to B. I. Shklovskii for the inspiration for this work and for numerous 
helpful suggestions. This work is supported by NSF under Grant DMR-9321417. 
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FIGURES 

FIG. 1. The quasiclassical image of the bubble phase, (a) Top view. The bubbles (dark 

circles) arc the places where the accumulation of the guiding centers occurs, (b) The enlarged view 
of one bubble. The dark region shows the guiding center density h'{x,y)/2TrP, while the toroidal 
figure illustrates the charge density distribution pj^{x,y) around the bubble. Half of the charge 
density is removed. This charge density is created by electrons moving in the cyclotron orbits 
centered inside the bubble. One of such cyclotron orbits is shown by the arc with the arrow. 

FIG. 2. The cohesive energy of the CDW as a function of for different numbers of electrons 
in a bubble M. The calculations are made for A?^ = 5 and = \/2. The crosses show the Laughlin 
liquid energies. 
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TABLES 



z^Af = 1/3 

1 -0.0853 -0.0819 -0.0733 -11.7% 

1 1 -0.1692 -0.1824 -0.1726 -5.6% 

2 2 -0.1970 -0.1939 -0.2163 10.3% 

3 3 -0.2135 -0.2056 -0.2433 15.5% 

4 4 -0.2251 -0.2228 -0.2480 15.1% 

5 5 -0.2341 -0.2367 -0.2767 14.5% 

VN = 1/5 

N M EU^L E^,, E^r ^E/E^r 

1 -0.0396 -0.0639 -0.0622 -2.7% 

1 1 -0.0986 -0.2109 -0.2043 -3.2% 

2 1 -0.1164 -0.2465 -0.2454 -0.5% 

3 2 -0.1267 -0.2474 -0.2811 12.0% 

4 2 -0.1340 -0.2538 -0.2990 15.1% 

5 3 -0.1395 -0.2562 -0.3187 19.6% 

TABLE I. The cohesive energies of the Laughlin hquid E^^^ and the CDW E^^^ in the 
Umit Nvg — ^ for different N. M is the optimum number of electrons per bubble. All 
the energies are given in the units of rghuc- The right column shows their relative difference: 
(-^coh ~ -^coh^) /-^coh^- "^^^ energy per electron in the uniform uncorrelated state E^^^ is pro- 
vided for reference. 
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VN = 1/3 

1 -0.1206 -0.1158 -0.1037 -11.7% 

1 1 -0.1297 -0.1520 -0.1424 -6.7% 

2 2 -0.1136 -0.1141 -0.1188 4.0% 

3 3 -0.1034 -0.0948 -0.1018 6.9% 

4 4 -0.0965 -0.0822 -0.0896 8.2% 

5 5 -0.0914 -0.0743 -0.0805 7.7% 

UN = 1/5 

N M E^^^ E™ ^E/Egpw 

1 -0.0560 -0.0903 -0.0880 -2.7% 

1 1 -0.0765 -0.1731 -0.1692 -2.3% 

2 1 -0.0677 -0.1423 -0.1396 -1.9% 

3 2 -0.0618 -0.1143 -0.1202 4.9% 

4 2 -0.0577 -0.0979 -0.1050 6.8% 

5 3 -0.0547 -0.0856 -0.0946 9.7% 

TABLE II. Same as Table | but = \/2, which corresponds to the electron density of 
1.6 • lO^^cm^^ in GaAs-GaAlAs heterostructures. The energies are now given in the units of 
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FIG.l. Fogler and Koulakov, PRL, "The absence..." 




FIG.2. Fogler and Koulakov, PRL, "The absence..." 



